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Abstract 

We discuss a new method of integration over matrix variables based on 
a suitable gauge choice in which the angular variables decouple from the 
eigenvalues at least for a class of two-matrix models. The calculation of cor- 
relation functions involving angular variables is simple in this gauge. Where 
the method is applicable it also gives an extremely simple proof of the clas- 
sical integration formula used to reduce multi-matrix models to an integral 
over the eigenvalues. 
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Matrix models have stimulated a great deal of interest in recent years following 
the discovery |l] that they provide a powerful tool in studying two dimensional grav- 
ity and its coupling with c < 1 conformal matter beyond perturbation expansion 
in the genus 0. In spite of that there has been very small progress in the actual 
mathematical tools used in solving matrix and multi-matrix models. In particular 
the multi-matrix models for which it is known how to perform the integration over 
the angular variables are still the ones listed in ref. ||. In all these cases the 
integration over the angular variables is done by using one fundamental formula: 

l f -r\ {1) \ {2) \ 

n(n-i)/2 i det e CAa A b 
M } \c) N\ A(A«)A(A( 2 )) [) 

where U is a unitary N x N matrix, A^ 1 -* and A 1 - 2 -* two diagonal matrices (A^ = 
<UA«)and 

A(A (J) )= II (Ai J) -A?) (2) 

l<b<a<N 

In eq. (|l|) d^U) denotes the Haar measure of SU(N) with the following normal- 
ization: 

This equation was first derived, for a general compact Lie group, by Harish- Chandra 
HJ , rediscovered in the context of matrix models in || and fully exploited in j3|. 
For a general proof of (JTJ) see also 0. 

In this paper we will present an alternative way to perform the integration over 
the angular variables, which consists in fixing the gauge in such a way that the 
integration over the remaining variables is trivial. When it is applicable (we shall 
discuss this point later) this method offers the advantage that after the gauge fixing 
the angular variables are decoupled from the eigenvalues, so that the calculation 
of correlation functions involving angular variables is greatly simplified. Also when 
such decoupling does not occur, as in the quartic couplings of the six vertex model 
(see page 26 in Ginsparg review paper 0), the dependence of the action on the 
angular variables is much simplified and it leaves some hope that the integration 
over them might eventually be done. Finally the method provides (in the cases 
discussed below) an extremely simple algebraic proof of eq. ([]]) . 

As a preliminary example let us consider the case of a one matrix model defined 
in terms of a hermitian N x N matrix A by the partition function 

Z N 



dm = K Z N _, d = f n N (3) 



dAe -Tr[V(A)] (4) 

The action TrV(A) is invariant under the unitary transformation 

A -> A' = TJ- l AU (5) 
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and so it depends only on the eigenvalues A a of A. We can treat the invariance 
under (||) as an ordinary gauge invariance and eliminate the unphysical degrees of 
freedom by a gauge fixing procedure Q In the gauge where A is diagonal 

K = 5 a b \ a (6) 

an infinitesimal unitary transformation of parameters e b is given by 

5A a b = z[e,A} a b = e a b (\ a -\ b ) (7) 

and the corresponding Faddeev-Popov determinant is 

det-^ = \{{\ a - \ b f (8) 

d€ c a>b 

The partition function Z N is then reduced to an integral over the eigenvalues: 

Z N = Q N f II d\ a II ( A « - ^fe- y(Aa) (9) 

J a a >b 

where Qn is the volume of the gauge group SU(N) defined in (|3|). 

Let us consider next the case of a two matrix model defined by the partition 
function: _ 

Z% ] = [ d V d^e- Tr ^ +v ^ +c ^ (10) 



where (p and <pft are respectively a complex N x N matrix and its hermitian 
conjugate. If one defines Lp = A + iB this is clearly equivalent to a model with 
two hermitian matrices A and B. The action in ([H]) is invariant under unitary 
transformations of ip: 

(p-xp = tftpU , <p j -> </ = £/y £7 (11) 



It is always possible , by using the invariance (|TTD , to reduce (p to a triangular 
form, that is: 

ip ab = a > b , <pl b = a < b (12) 

Notice that in such gauge the diagonal elements of ip and ^ coincide with the 
respective eigenvalues: ip aa = X a and (p\ a = X a . The advantage of the gauge fixing 
conditions fll2"|) is that the corresponding Faddeev-Popov determinant is very simple 
and depends only on the eigenvalues on <p and ip* . In fact if one denotes by 
the l.h.s. of the gauge conditions (|12|), that is Q ab = p ab for a > b and § ab = <p\ b 
for a < b, and by e ab the parameters of an infinitesimal unitary transformation 
(e&a = Sob), then it is easy to show that: 

det 5 -pL = A(A)A(A) = J] |A a - A fe | 2 (13) 



8e cd 



a>b 



: A number of people must have been aware of this procedure, to my knowledge it has appeared 
in the literature in Ginsparg review paper and in H . 
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With this choice of gauge the partition function (JTD) takes the form: 

f N _ 

4^ = V N J n d\ a d\ a II d^abdWab A(A)A(A) 
a=l a<6 

iV 

ex P {- Et^(Aa) + V(\) + CAX] - C^abWab) (14) 
a=l a<fe 

As the angular variables appear quadratically in ([14]) and are decoupled from the 
eigenvalues , the corresponding gaussian integral is trivial and gives: 

Z%> = (—) Q N / II d\ a dT a A(\)A(\)e-^lvM+vix a)+ cX a A a ] (15) 

^ C ' J a=l 

Eq.(^) is the same as the the one considered for instance in || , except that the two 
hermitian matrices are replaced here by complex conjugate matrices ip and <p>. The 
novelty of our procedure is in the simple quadratic dependence of the action from the 
remaining angular variables after fixing the gauge as shown in eq. (0) . This allows 
in particular to calculate (or at least to reduce to integrals over the eigenvalues) 
correlation functions of gauge invariant quantities , such as for instance Trip 2 ^ 2 
and TrLpip<(p(p\ that involve also angular variables. In these cases the integration 
over ifab and Tp ab can be done by simple Feynman diagrams techniques. These 
correlation functions are the first terms of the expansion in powers of the coupling 
constants of the partition function generated by the action 

S = Tr[V((p) + V(ip^) + api^ + g\<p<p<p*$ + ^wVv^ (16) 

which contains as a particular case the matrix model describing the coupling of the 
six vertex model to gravity. In ( |T5D the angular variables are not decoupled , in 
fact in the gauge (|i~2|) we have for instance: 



Tripip^ipip^ = ^ ||A a | 4 + 2|A a | 2 \(pba\ 2 + \¥ab\ 2 + 2\ a (p ba ip ac ip bc 

+2\ a tp ab tp ca !p cb + PabPcdVcbPad} i 17 ) 

where the Safe's and the Tp ab s at the r.h.s. of the equation stand for the off-diagonal 
elements only ( ip a b = Tp ab = for a > b). In this situation the task of reducing the 
calculation of the partition function to an integral over the eigenvalues is , beyond 
the first few terms in the coupling constants expansion, a formidable one but not , 
according to some preliminary calculations, a completely hopeless one. 

It would obviously be very interesting to generalize these results to the case of 
models of two hermitian matrices. However a hermitian matrix cannot be cast in a 
triangular form and we could not find any other gauge choice that offers the same 
advantages. Still one can apply our method to a different model, which seems to 
be closely related to the one of two hermitian matrices although we were unable to 
prove a complete equivalence. Let us consider two hermitian matrices Ai (i = 1,2) 
parameterized as usual in the following way: 

Ai = UiKU} (18) 
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where Aj are diagonal (Aj ab = Sg^X^') and Ui are unitary matrices. Let us associate 
to any hermitian matrix (|i"8|) a real matrix A\ denned as follows: 

A = a^A 1 (19) 

where Ui is the matrix of SL(N,r) obtained from Ui by a "Wick rotation" of the 
group parameters corresponding to the real- symmetric generators of SU(N). More 
explicitly, if we denote by the antisymmetric generators of SU (N) and by 
the symmetric ones, we have: 

JJ, = e *r$T(<*)+UT$S°* °ab^i*a b = e *r<$T(<*)-*$S<* ^0) 

It should be noticed that the A^'s are not generic real matrices since their eigen- 
values are constrained to be real. By means of a SL(N, r) transformation the two 
matrices A^ can be cast into a triangular form, more precisely: 

A2=0 a>b ii?=0 a<b (21) 

As before the diagonal elements coincide in this gauge with the eigenvalues. Con- 
sider now the partition function 

Z { n ] = I d^dJ^e~^ m ^ m ^ wm \ (22) 



where the integration is restricted to a domain T> corresponding to real eigenvalues 
for the real matrices A^> and A^> . The action is invariant under transformations 



of the non compact group SL(N, r) ; by fixing the gauge according to ([H]) and by 
calculating the Faddeev-Popov determinant in complete analogy to eq. flT5| ) one 
finds: 

, N 

Zf = Cl N J dX^dX? n dA^dAf} A(A«)A(A( 2 >) • 

a=l a<b 

exp(-E [Vim + V 2 (\W)+ic\M\W} -ic^AWAw) (23) 

[ 0=1 a<b J 

where Qn is the infinite volume of the gauge group. Notice that the coupling 
constant c in front of the term TrA^'A^- 2 ' is multiplied by i in (p^). Such "Wick 
rotation" is needed in order to have a finite result from the integration over the 
remaining angular variables in ( |2"3"D which would otherwise be divergent due to the 
non compact nature of the group. The final result is: 

/9 „ N N(N-l)/2 r N 

4 2) = CIn(^) /n^^Ai 2 )A(A«)A(A( 2 )). 

V C J J a=l 



e -E. 



(24) 
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This result coincide (after substituting ic back with c) with the standard result 
obtained in || for hermitian matrices. However it should be noticed that , as in 
the case of complex conjugate matrices , eq. ([23]) gives some extra information, 
namely that the angular variables are decoupled from the eigenvalues and appear 
quadratically in the action. All the considerations following eq. fllil) and (|15|) can 
be applied to the present case by simply replacing ip and ip^ with A^ 1 ' and A^ 2 \ in 
particular eqs. similar to (|T6| ) and fllTD can be easily derived. 

It is possible also to prove the SL(N,r) analogue of eq. ([[]), which is in fact a 
simple consequence of (|23D and of 



dA (l) = J] d\f& 2 {\ {i) )dii{Ui) (25) 

a 

By inserting ([19]) and fl25|) into eq. (^2|) one finds: 

Z N = Cl N [Y[dX^d\^d^U)A\\^)A\\^). 

e -E a [vi(Ai 1) )+^(Ai 2) )]- l cTr(t/A(i)C/-iA(2)) 



It is now sufficient to compare eqs (|26|) and (^3|) and to take into account that the 
integral over the eigenvalues is weighted by arbitrary functions which are symmetric 
under separate permutations of the eigenvalues of A^ and A^ 2 \ The result is: 

N\ J d»(U)e- icTr & Am - lA ^ = 



V A(A«)A(A( 2 )) J ^ ab ba 
N{N-i)/2 det I e JCAa A & J 

"7 J a(a«)A(a( 2 )) (27) 

where the sum is over all permutations P on the eigenvalues A^. In deriving 
eq.(p7[) one has used the fact that the l.h.s. is invariant under permutations of the 
eigenvalues. In fact it is easy to show that for any permutation P it exists a matrix 
Up G SL(N, r) such that 

PA« = UpA^Up 1 (28) 



Eq. fl27|) can be generalized by including in the original partition function ( P2|) 
higher order interactions involving A^ and A^ 2 \ For instance the addition of a 
quartic term igTrA^A^A^A^ to the action in (^) would result in an extra 
term 

- ig (TrUA^U^A^UA^U- 1 ^) (29) 

in the exponential at the l.h.s. of (|2Tf) . Correspondingly the exponential at the 
l.h.s. would acquire a contribution of the type 

-i9 {E [Ai 1} a A?> 2 + 2A«A( 2 ) (A«i£) + + 2A«ii?i«i( 2 a ) + 

+2 a( 2 U {1) 4 (2) 4« + a {1) a (2) a {1) A {2) ]\ rsm 

+ Z A a ^ab A bc ^ca + ^afe ^fec ^ccZ ^da J J l JU J 
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where A® denote the off-diagonal elements only and the sum is then restricted to 
Aab = f° r a > b an d ^ab = f° r o, < b. Clearly in this case the integration over 
A^b and A^} at the r.h.s. of eq. (27) involves the eigenvalues and it cannot be done 
explicitly, at least so far. 

It should be noticed that our proof of eq.([27]) and its generalizations follows in a 
purely algebraic way from comparing two different choices of gauge in the partition 
function. 
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